Optimal probabilistic cloning and purification of quantum states 
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We investigate the probabilistic cloning and purification of quantum states. The performance of 
these probabilistic operations is quantified by the average fidelity between the ideal and actual output 
states. We provide a simple formula for the maximal achievable average fidelity and we explictly 
show how to construct a probabilistic operation that achieves this fidelity. We illustrate our method 
on several examples such as the phase covariant cloning of qubits, cloning of coherent states, and 
purification of qubits transmitted via depolarizing channel and amplitude damping channel. Our 
examples reveal that the probabilistic doner may yield higher fidelity than the best deterministic 
doner even when the states that should be cloned are linearly dependent and are drawn from a 
continuous set. 
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I. INTRODUCTION 

The recent spectacular development of the quantum 
information theory has revealed that information pro- 
cessing based on the laws of quantum mechanics enables 
implementation of tasks that are impossible or very hard 
to accomplish classically. The prime examples are the un- 
conditionally secure cryptography p| and the exponential 
speedup of certain computational tasks, such as the fac- 
toring of integers On the other hand, the linearity 
of quantum mechanics also imposes certain constraints 
on the processing of quantum information that have no 
classical counterpart. Perhaps the most famous example 
is the no-cloning theorem which states that an unknown 
quantum state cannot be copied Q. However, this re- 
striction provides, in fact, a valuable resource explored 
in the quantum key distribution protocols, because it 
forbids an eavesdropper Eve to gain information on the 
distributed secret key without introducing errors. 

Since exact copying is forbidden, a natural problem 
arises what is the optimal approximate cloning trans- 
formation. This question was first asked by Buzek and 
Hillery in their seminal paper Q and since then it has 
been addressed by numerous authors who considered var- 
ious cloning scenarios, such as cloning of qubits [a. |g, 
cloning of d-dimensional systems (qudits) lEM [lO, [ij 

and cloning of continuous variables [niinij Much 

attention has been recently paid also to the cloning of 
subsets of Hilbert space, such as the phase covariant 
cloning machine for equatorial qubits [1^ [la, [l3> Ull and 
cloning of real states [T^ or maximally entangled states 
[23 | . Remarkably, the cloning machines turned out to be 
very efficient or even optimal eavesdrop ping attacks on 
many quantum cryptographic protocols |2lL 123. 123. l2^ 
and their investigation is largely motivated by these prac- 
tical aspects. 

Typically, the doner is assumed to be a deterministic 
machine that always produces an output. Nevertheless, 
one can consider also probabilistic cloning machines that 
sometimes fail and do not provide any outcome. The 
probabilistic doners have been discussed in the litera- 



ture in the context of cloning of a discrete finite set of 
quantum states and it was shown that an exact proba- 
bilistic cloning is possible if and only if the set consists 
of linearly independent states [2^ |2g . 

However, one may hope that the probabilistic machines 
may yield better results also for the sets of linearly de- 
pendent quantum states and even for infinite (continu- 
ous) sets. Here, we investigate the probabilistic cloning 
of linearly dependent states and we establish a general 
theory of the optimal probabilistic cloning machines. We 
provide a simple formula for the optimal average fidelity 
of the probabilistic machine and we also show how to 
construct the optimal cloning transformations. 

In fact, our formalism is very general and it con- 
cerns optimal probabilistic implementations of arbitrary 
transformations whose outputs should ideally be pure 
states. Besides cloning, this includes also universal NOT 
gate for qubits l28l|. and, perhaps even more impor- 
tantly probabilistic purification of mixed quantum states 
|29ll30l l3ll | . In what follows we first establish the general 
formalism and then we work out several explicit examples 
that will illustrate our method. 



II. OPTIMAL PROBABILISTIC 
TRANSFORMATIONS 

The probabilistic machines investigated in the present 
paper optimally (in a sense defined below) approximate 
the map from a set of input (generally mixed) states 
Pin to the set S'out of the output pure states iV'out), 



Pn 



^out(Pii 



(1) 



where i/'out = |''/'out)('0out| is a short hand notation for 
the density matrix of a pure state. If S'in is a set of pure 
states (such as in the case of cloning) then we replace 
with TpiT^. The most general probabilistic transformation 
in quantum mechanics is a linear trace decreasing com- 
pletely positive (CP) map [s^l that transforms operators 
on the input Hilbert space Tiin onto the operators on 
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the output Hilbert space Hout- According to the Jami- 
olkowski isomorphism [s^, any CP map pout = ^(Pin) 
can be represented by a positive semidefinite operator E 
on the Hilbert space H = Hin "^out- Let |j) denote a 
basis in a d-dimensional Hilbert space Tijn . The operator 
E can be obtained from the maximally entangled state 
on Hin (8) Ti-in, \^) = J2j=i \j)A\j)B, if we apply the map 
£ to one part of |$). We have E = Ia® Ssi^) where I 
denotes the identity map. 

The transformation pout = ^(Pin) can be rewritten in 
terms of E as follows, 



Pout 



TiiniEpl'^tc 



(2) 



where Trin denotes the partial trace over the input 
Hilbert space and T stands for the transposition in the 
basis The map must be trace decreasing which means 
that Tr[pout] < Tr(pin) for all pin > 0. This implies that 
E must satisfy the inequality 



Trout < li: 



(3) 



where 1 denotes the identity operator and the equality 
in (PI is achieved by deterministic (trace preserving) CP 
maps. 

Consider a particular input pm- The normalized out- 
put state Pout is given by pout = Pout/-P(pin) where 

P(pi„) = TiiEpl ® lout] (4) 



is the probability of successful application of the map £ 
to Pill. The performance of the map £ for the particular 
input Pill can be conveniently quantified by the fidelity 
between the actual and the ideal outputs: 



F(pin) = (l/'outlPoutlV'out) 

Expressed in terms of E we have 



Fipin) 



1 



TviEpl^(g,ib,nt). 



(5) 



(6) 



We assume that the set S-m is endowed with an a-priori 
probability distribution dpin such that Jg, dpin = 1. 
Here and in what follows we assume that the set Sin is 
continuous. Of course, all formulas remain valid also for 
discrete sets, one simply has to replace the integrals with 
corresponding summations over the elements of Sm- 

The average probability of success is defined as 



where 



P= / P(piii)dpin = Tr[SA], 



A= I (g) lout dpiii. 



(7) 



(8) 



We now introduce the mean fidelity F of the transforma- 
tion £ as the average of the fidelities F(piii), with proper 
weights P(pin)dpiii/P, 



F : 



F[Pin) — ^—apu 



(9) 



The mean fidelity is the figure of merit considered in the 
present paper and in what follows we shall look for the 
optimal map £ that maximizes F. 

If we insert the expression (j^J) into Eq. © we obtain 
F = F/P where F = Tr[PP] and 



^= Pin i^outdpu 

"'Sin 



(10) 



Taking everything together, we want to find E that max- 
imizes the mean fidelity 



F 



Tt{ER) 
Tt{EA) ' 



(11) 



where P > and A > are defined above. The positive 
semidefinite operator E representing a trace-decreasing 
CP map must satisfy the constraint (PJ. However, this 
constraint is irrelevant as far as the mean fidelity Hll|l is 
concerned because the value of F does not change under 
the re-normalization 



E 



E' = e^^ E 



(12) 



where Cmax — max[eig(Trout-E')] and E' satisfies the in- 
equality © by construction. This fact greatly simplifies 
the analysis. Strictly speaking, these arguments are valid 
only for finite dimensional Hilbert spaces where Cmax is 
always finite and P' = Tr[P'^] > since A > 0. As we 
will see in the next section, a little extra care is needed 
when dealing with infinite dimensional systems. 

The above argumentation shows that we have to max- 
imize the fidelity (|ll|l under the constraint £^ > 0. With- 
out loss of generality, we can assume that the optimal 
P is a pure state E — |P)(P|. We introduce new state 
IP) = A^I'^\E) and rewrite as follows. 



P = 



{E\A-^I^RA-^I^\E) 
{E\E) 



(13) 



It follows that the optimal vector \E) is the eigenvector 
l^max) of M = A"i/2pA"i/2 that corresponds to the 
maximal eigenvalue /Xmax of M . The maximal achievable 
mean fidelity is equal to the maximal eigenvalue, 

Pinax = max[eig(M)] = max[eig(A-ip)] . (14) 

This formula is one of the the main results of the present 
paper. The transformation that achieves Pmax is ex- 
plicitely given by 



^ — ^max^ ^^^|Mmax)(Mmax|^ ^^'^ , 



(15) 



Sin 



where we have normalized according to 1)12(1 so that E 
is a trace-decreasing map. If the largest eigenvalue /Zmax 
is non-degenerate, then this is the unique optimal E and 
the problem is thus completely solved. However, if the 
eigenvalue Pmax is n-fold degenerate, with |pmaxj)j j = 
1, . . . , n being the n eigenvectors, then there exist many 
different transformations that saturate the fidelity bound 
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(|14(l . It can be proved by direct substitution into Eq. 1111) 
that any operator 

n 

E,kA-^'^\^Jin..y,.,J){^ira.^,k\A~^'^, (16) 

3,k=l 

yields the maximal fidelity F = /imax- Let K, be the 
Hilbert space spanned by the vectors A~^/^|yUinaxj)- 
Then E can be any positive semidefinite operator on /C 
that satisfies Q. In this case, we would like to find the 
map E that maximizes the average probability of success 
P while reaching the fidelity i^max- The optimization 
problem that has to be solved can be formulated as fol- 
lows, 

maximize P — Tv[EA\ under the constraints 

(17) 

where B{JC) denotes the set of linear bounded operators 
on /C. This is an instance of the so-called semidefinite 
program (SDP) that can be very efficiently solved nu- 
merically and by means of the duality lemma one can 
easily check that the global maximum was found 'S^j. 
In this context it is worth noting that many optimiza- 
tion problems in quantum information theory can be for- 
mulated as semidefinite programs. This includes several 
separability criteria 35, 36], calculation of distillable en- 
tanglement fSVl f38l | . determination of optimal POVM for 
discrimination of quantum states [s^ llOj] , derivation of 
optimal trace-preserving CP maps for clonin g Ell |4^ . 
construction of local hidden variable theories [43 etc. 

Generally, the fidelity -F(pin) will depend on pi^- How- 
ever, there is an important class of sets of input states Sin 
and transformations ^ such that the optimal CP map is 
universal. By universality we mean that the probability 
of success P{pin) as well as the fidelity F{pin) is indepen- 
dent of Pin . This occurs whenever the set of the input and 
output states can be obtained as orbits of some group G. 
Consider a compact group G with elements g. Let U{g) 
and V{g) denote unitary representations of G on Hin and 
Hout, respectively. The unitary U{g) generates the set of 
input states, 

Pinig) = U{g)pMUHg), (18) 

where go is the identity element of the group and U (go) = 
1. We also assume that the set of output states can be 
obtained from '!/'out(<?o) as follows, 

^out(<?) = ^(5)^out (50)^^(5) (19) 

and V{gQ) — 1. The final assumption is that the distri- 
bution of the inputs coincides with the invariant measure 
on the group G, dp\n — dg. Under these assumptions, it 
is possible to convert any optimal map £ into a universal 
map £ that achieves the same fidelity as £ by the twirling 
operation. One first applies randomly a unitary U{h) to 
the input and then this is undone by applying V~^{h) to 



the output. The composition of the twirling operation 
with the map £ yields 

Pout(5)= / V\h)i:r,n[Epl,{hg)®t,nt]V{h)dh (20) 

JG 

and the probability of success reads 

P'[pU9)] - / "^^[EpLihg] ® ^ont]dh = P. (21) 

JG 

Here, we used the group composition law U{h)U{g) — 
U (hg) to obtain U{h)pin{g)U^h) = pin{hg), and the sub- 
stitution q — hg, dq = dh. Similarly, we find that 

^'(Pin) = i / ^4Epl{hg) (g> ^outihg)]dh = F, (22) 

^ JG 

which confirms that the twirling operation results in a 
universal machine that works equally well for all possible 
input states. 

III. PROBABILISTIC CLONING 

Having established the general formalism, we now turn 
our attention to the explicit examples of application. Let 
us first consider the universal symmetric 1 ^ M cloning 
machine for qubits. Here, the input state is a single qubit 
= cos^|0) + e*"^sin||l), uniformly distributed over 
the surface of the Bloch sphere, d-f/'in = sinddddtj)- 
The cloning machine should produce M identical clones, 
hence V'out = "^in*^ ■ The operators R and A can be easily 
calculated with the use of the Schur lemma, 

R ^ j^{Tl+M+if\ A=Ui®n+,M, (23) 

where H+^m denotes a projector onto symmetric sub- 
space of M qubits and Ti indicates partial transposition 
with respect to the first qubit. On inserting the opera- 
tors ^ into Eq. 13 we find that i^max = 2/{M + 1). 
As shown in Ref. |7| the optimal deterministic cloning 
machine saturates this bound, hence it is impossible to 
improve the fidelity via probabilistic cloning. 

Let us now consider the transposition operation for 
qudits, i.e. a map that produces a transposed qudit state 
ijj'^ (in some fixed basis) from N copies of ip, -0*^^ 
ijj^ . For qubits, this map is unit arily equivalent to the 
universal NOT gate ^ V'-L iS El ■ The Hilbert 

space Tiin is the fully symmetric subspace of the Hilbert 
space of qudits and Tiout is the Hilbert space of a 
single qudit. In the formulas (jHl and Hl()|l for A and R 
we average over all V'in that are represented as orbits of 
the group SU{d) according to Eq. ((TH)l . The probabifity 
density dipin = dg^ where dg is the invariant measure on 
the group SU[d). With the use of the Schur lemma one 
easily finds 
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where n^''jY^j^ is the projector onto symmetric subspace 

of TV + 1 qudits and D{N + l,d) = is the dimen- 

sion of this subspace. The optimal fidelity obtained from 
Eq. lO reads F^ax ^ {N + 1)/ {N + d), which is exactly 
the fidelity of the optimal deterministic estimation of the 
qudit state from N copies 44|. Note that the fidelity 
-Fmax = 2/((i + 1) of the optimal deterministic transposi- 
tion map for N = 1 was recently derived in |45l | . 

In all the above examples the optimal probabilistic ma- 
chine could not outperform the deterministic machines. 
This can be attributed to the very high symmetry present 
in all the above considered examples. The question is 
whether there are interesting cases when S-^ is a contin- 
uous set of linearly dependent states and the probabilis- 
tic machine achieves higher fidelity than the deterministic 
one. Below we answer this question in affirmative by pro- 
viding explicit examples. We will focus on phase covari- 
ant cloning machines, where the underlying group is the 
Abelian group U{1). Specifically, we shall first consider 
probabilistic N ^ M phase covariant cloning of qubits 
[11 [13. Here, the input state ji/') = (|0) + e*'^|I))/%/2 
lies on the equator of the Bloch sphere and is charac- 
terized by a single parameter, the phase (j>. Moreover, 
J dg = J^^ d<l) / {2t:) . The input and output states are 
given by 



(25) 



The integrals appearing in the expression (|I()|I can easily 
be carried out and one arrives at 



1 ^ 

^ " 2M+N \^M,N,y){^M,N,yl 



(26) 



V=-N 



where 



min(N,M~y) 



1$ 



M,N,yl 



E 



fc— max(0,— y) 




\N,k)\M,y + k), 
(27) 

and I A^, k) denotes a totally symmetric state of N qubits 
with k qubits in the state |1) and N — k qubits in the 
state |0). Similarly, one gets 

I ^ fN\ 

^ = ^Eu)'^'^^^^'^'^^°^'- ^^^^ 

fc=0 



Since the states \^M,N,y) are mutually orthogonal, the 
matrix R is diagonal and, consequently, also M = 
A~^^'^RA~^/'^ is diagonal. The maximal eigenvalue can 
thus be easily determined and the maximal fidelity is 
given by 



1 ^ 

k=0 



M 



(29) 



where [x] denotes the integer part of x. For > I, 
the fidelity (|29|l is higher than the fidelity of the optimal 



deterministic phase covariant doner that was given in 
[l7| . The improvement of the fidelity is typically of the 
order of one percent. 

The optimal probabilistic cloning transformation can 
be written as 



|7V, fc) 



M 



N 



Af]] \k + AmnJ V k 
where Amn = [^^] and 



J\f ~ max 
fc 




|M, AuN + k), 

(30) 



(31) 



is a normalization prefactor. If Af — is even, then Eq. 
(|30|l is the unique optimal phase-covariant probabilistic 
cloning transformation that optimally matches the input 
state IV-)®^ onto the ideal output For odd M - 

Af, however, we can obtain another optimal operation by 
replacing A jvfjv with Am AT + 1 in Eqs. ijSU)! and (|?T|l . This 
implies that the support K. of the optimal operator E is 
two dimensional and the optimal map that maximizes 
the success probability has to be calculated by solving 
the semidefinite program (|17|l . 

A second example where the probabilistic doner out- 
performs the deterministic one is the 1 ^ M copy- 
ing of coherent states \a) on a circle. Recall that 

|a) = e^l"l S^^o Q^^/"^/"^ I")' where \n) is the n- 
photon Fock state. Since we assume that \a) lie on a 
circle, the amplitude r = \a\ is fixed while the phase 
(j) — arg(a) is arbitrary. First of all, we observe that 
the perfect cloning is equivalent to noiseless amplifica- 
tion, because the output state ja)®*^ can be unitarily 
mapped onto the state \VMa) (g) |0)®*^~^ by an array of 
Af — 1 beam splitters . Thus the cloning is equivalent 
to I a) |\/Afa). The operators A and R are calculated 
as averages over the phase (j) and one finds that the max- 
imum eigenvalue of Af is /imax — 1 which indicates that 
an exact probabilistic cloning is possible. However, we 
must be careful because we deal with infinite dimensional 
Hilbert space and it turns out the the fidelity F = 1 can 
be achieved only in the limit of zero probability of suc- 
cess P ^ 0. Nevertheless, arbitrarily high fidelity can be 
reached with finite success probability if we first project 
the input state onto the subspace spanned by the first 
A^ -I- 1 Fock states |0), . . . , jA^) and then apply a diago- 
nal filter that approximates the noiseless amplification, 
|n) ^ Af("-^)/2|n), n = 0,---,Af. The fidelity 



N 



F 



Af"|a|' 



ri=0 



,1 



(32) 



can be arbitrarily close to 1 as A^ ^ oo. Clearly, this 
probabilistic cloning achieves even higher fidelity for the 
coherent states inside the circle, where \a\ < r. 
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IV. PURIFICATION OF MIXED STATES 

Another important application of the optimization 
technique developed in Sec. II consists in the design of 
the optimal protocols for purification of mixed quantum 
states. Suppose that Alice and Bob can communicate via 
a noisy quantum channel C. Alice wants to send to Bob 
a quantum state tp from some set Sin- However, since the 
channel is noisy, Bob receives mixed state p — C{tp). To 
partially compensate for the effects of the noisy channel, 
Alice sends N copies of the state ip to Bob who subse- 
quently attempts to extract from the state p®^. The 
purification of qubits transmitted through the depolariz- 
ing channel, p = ml' + ^(1 — ?7)1, has been analyzed in 
detail in Refs. [23, H^]- Very recently, the optimal pu- 
rification protocol for two copies of the qubit has been 
demonstrated experimentally for the polarization states 
of single photons by exploiting the interference of two 
photons on a balanced beam splitter [^J- Applications 
of the purification procedure to the quantum state es- 
timation and transmission have been discussed in Refs. 
[46l f47| . Here, we demonstrate that the present opti- 
mization procedure can be used to straightforwardly de- 
termine the optimal probabilistic purification protocol. 
Then we will consider the same problem for the ampli- 
tude damping channel. 



A. Depolarizing channel 

For the sake of simplicity we illustrate the method on 
the case when Alice sends two copies of the state 1-0) to 
Bob, hence Bob's input mixed state reads pin = {rjip + 
i(l — 77)1)®^. The ideal output is a single-qubit pure 
state ip. Assuming uniform distribution of ip over the 
surface of the Bloch sphere, one obtains the following 
expressions for the operators A and R, 



0.06 



A 



^H+^123 



1 - 77^ 



L123 



1, 



- 7/)(ai ® n+,23 + 12 ® n+,13), (33) 

6 

where 1 and 2 label the input qubits while 3 labels the 
output qubit, Tl+.jk and H+ ij^ are projectors on the 
symmetric subspace of two qubits i,j or three qubits 
i,j,k, respectively, and Ta stands for the partial trans- 
position with respect to the third qubit. From Eq. 114|) 
we obtain the maximal achievable purification fidelity: 



F 



3 + 477 + 772 
2(3 + 772) 



(34) 



which is larger than the original fidelity Fq = (^\MA_ = 
(l+77)/2,forallO < 77 < 1. As shown in Refs. [29l.l3al3ll|. 

the optimal purification strategy is to project the two- 
qubit state onto the symmetric subspace and then 



0.04 



< 

LL 

< 



0.02 



0.00 




FIG. 1: The maximal possible improvements of fidelity AF{ri) 
(solid line) and AF± (77) (dashed line) that can be achieved by 
purification when the two-qubit state IV'i/') or It/ji/jx), respec- 
tively, is sent through the depolarizing channel with parame- 
ter 77. 



throw away one of the qubits. This procedure achieves 
the optimal fidelity (|34|l . Let us now demonstrate that 
this protocol can be derived by solving the semidefinite 
program (|17|l . The maximal eigenvalue of the matrix M 
is doubly degenerate and the basis states that span the 
two dimensional Hilbert space K, are given by 

|ei) = i=(V2|vI/+)|0)+2|ll)|l)), 
|e2)--^(V2|vI/+)|l)+2|00)|0)), (35) 

where |*+> = (|01> + |10))/V2. _ From Eq. and 
the definition of A wc find that P = Trin[AX], where 
A = Js.^pfndPin and X = Trout (-E)- Since E = 
^j=i 2 Cjfe|ej)(efc| it follows that the support of X is the 
symmetric subspace of two qubits. From Eq. (PJ we thus 
have X < H_|_.i2. The maximum P is obtained when 
X = n+,12, which can be achieved by the following choice 



E : 



:(|ei)(ei| + |e2)(e2|). 



It is easy to check that this trace-decreasing CP map 
indeed describes the projection of two qubits onto the 
symmetric subspace followed by tracing over the second 
qubit. 

The strategy to send the state \'4'i') is not the only pos- 
sible option how Alice can encode the state into the 
two qubits that she sends to Bob via depolarizing chan- 
nel. For instance, she can send him the state |-i/))|7/) 1 ), 
where (i/'IV'-l) = 0- As shown by Gisin and Popescu |27l| . 
the state ji/') can be estimated with higher fidelity from a 
single copy of the state than from a single copy of 

1 7/^-0). We cannot therefore a- priori rule out that sending 
the state iV'V'i) can be advantageous also in the present 
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context. If Alice sends then Bob's input mixed 

state reads 



Pill 



The calculation of the optimal purification fidelity is com- 
pletely similar to the case of sending IV'V')- The integrals 
^ and H10|l yielding the relevant operators A± and R± 
can be easily evaluated with the help of the substitution 

tjj^ = ! — ■(/' and we obtain 



13 



l + rj2 1 „ 

^-^ii 12 - -rn+,12 

— n+_i23 H ^ ll23 

+?(! + '7)12 ® n+,13 - J(i - 77)ii (g) n+,23, 

6 6 

(36) 

The maximal fidelity F±^ is determined as the maximum 
eigenvalue of the operator A~^Ri_. 

The results of numerical calculations are given in Fig. 
1. For comparison we plot on this figure the gains in 
fidelity AFj_ = F^irf) - Fo{t]) and AF = F{r]) - Fo{r]) 
achieved when Alice sends the state \ip'ip±) or lip'ip), re- 
spectively. We can see that as far as the purification 
is concerned, it is strictly better for any < 77 < 1 to 
send the state \ipip} than \'ip'ijj±). The non-zero values 
of AFl {•!]) clearly show that purification is possible also 
when Alice sends the state IV'V'x) but the fidelity im- 
provement is much smaller than when sending the state 



B. Amplitude damping channel 

To further illustrate the utility and universality of our 
optimization method, let us now consider a different class 
of noisy channels, namely, an amplitude-damping chan- 
nel that maps a pure state onto a mixed state 



Pad{^,<P) = 



2 2 i9 

ry cos^ ^ 
§sin?9e** 1 



■| sin-iJe 



9 2 i9 
if COS f 



(37) 



This channel may arise, for instance, when the qubit is 
represented by the ground and excited atomic states \g) 
and |e) where |e) can decay to \g) via spontaneous emis- 
sion. In order to preserve the covariance ifTI^ that guar- 
antees the universality of the optimal purification pro- 
tocol, we shall assume that Alice is sending to Bob TV 
copies of a state \tp{(j))) = (|0> + e"l'\l))/V2 that lies on 
the equator of the Bloch sphere, i.e. the set S'in consists 
of the states pI^sW = Pad (V2, </-). 

The operators Rad and Aad are obtained by integrat- 
ing over the phase 0, 



PBOsi 



Rad 



1 f^'" N 



-(/)) (g> VX^) d(l>. (38) 



We now prove that when determining the maximum 
achievable purification fidelity, we can assume that the 
optimal map £ is a composition of two maps, £ = SoVn, 
where the map Vjy projects the input state p'^qb onto 
the symmetric subspace of N qubits TI+^n and £ maps 
operators on H.+^n onto operators on the Hilbert space 
of single output qubit. By definition, the operators Aad 
and Rad commute with arbitrary permutation operator 
IIj that changes the order of the N input qubits. Sup- 
pose that E — |e)(e| is an optimal map yielding maximal 
fidelity Fmax which is equal to the maximum eigenvalue 

X/2 1/2 

of matrix Mad — A^^jj RadAj^j^ . As shown in Sec. II, 

the corresponding optimal eigenvector of M is related 

1/2 

to |e) as follows, = AJ^^\e). Let us now consider a 
symmetrized state \x) that we obtain from |e) by making 
a linear superposition of all N\ permutations of N input 
qubits. 



En. 



(39) 



1 /2 

It is easy to show that A^^jx) is an eigenstate of Mad 
with eigenvalue -Fmax- We have 



MadAK^ \x) 



= A'J^A^],RADJ2u,A-/^'\^,) 

j 

^ ^AD ^j^AD ^maxl^) 



(40) 



where we have used that both Aad and Rad commute 
with IIj. Thus the map X = |2;)(a;| achieves the maximal 
fidelity F^ax- It holds that Trout (^) e B{n+^N) which 
proves that we can restrict our attention to the maps £ 
when calculating the maximal fidelity of purified state. 

The calculations can be further considerably simpli- 
fied by the observation that the optimal map £ can be 
made phase-covariant 16, 17], that is, invariant under 
the twirling operation, 

'^u*{4>)®v{4>)EU^{(p)®v\4>), (41) 

zvr 



E = 



where t/((/))|7V, fc) = e*'=*| A^, /c), fc = 0, . . . , TV and V\j) = 
j = 0, 1. This imphes that the operator E can be 
expressed as a direct sum, 



N 



E^ ^ Ek, 



(42) 



fe=-i 



where the support of the operator Fj, is a two- 
dimensional Hilbert space Ti^, fc = 0, . . . , A^ — 1 spanned 
by |A,fc)|0) and \N,k + \)\l) and 



F_i = /i_i|A^,0)(Af,0|®|l)(l|, 
En+1 = pN+i\N,N){N,N\®\{)){Q\. 



(43) 
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The decomposition (|42|l implies that we can perform the 
optimization of each CP map E^. separately and then 
choose the k that yields the highest fidelity. It is easy to 
see that the trace decreasing CP maps E^i and E^+i 
lead to very low purification fidelity i so it is optimal for 
all N to set = ijln+i = 0. Without loss of gener- 
ality, we can assume that the optimal Ek is a rank-one 
operator, Ek = \Ek){Ek\, where 



\Ek) 



\N, k)\0)+aN,kiv)\N,k + l)\l) 



(44) 



The action of this operation can be understood as follows. 
First the A^-qubit state is projected onto two-dimensional 
subspace of H+^n spanned by |iV, k) and |A^, fc + 1) and 
then the following transformation is carried out. 



\N,k) ^ |0), 



\N,k + l) 



aN.k{v)\^) 



(45) 



The un-normalized density matrix of the purified qubit 
obtained by applying the map 145|) reads 



Pout 



"k,k '^N,k"k,k+1^ 
(N) id, I 12 (N) 

O'N.kCTlk'+ie '^ \aN,k\ CTfc+i./c-Hi 



(46) 



The relevant matrix elements 

a^T = {N,j\pf^in/2,0)\N,k) 
can be expressed in terms of a finite series. 



'k,k 



= 2"" 77 



(=0 



k-l 



'k,k+l — ^ 'I 
k 

X 



N2N-2k-l 



k + 1 



E 



N -k 

N-k ■ 
k + l-l 



maximal probability P of optimal purification. We have 
therefore carried out full calculations of the operators 
A AD and Rad for = 1,2,3 and determined the de- 
generacy of the maximal eigenvalue of matrix M . These 

1 /2 

calculations reveal that the optimal eigenvectors A^^je) 
of M all satisfy the relation Trout [e] £ B{T-L+^n) so we can 
in fact consider only the maps of the form 144|) without 
any loss of generality. 

If only a single qubit is sent to Bob, then the only 
option is fc = and the best filter obtained by setting 
Q^i,o = vl \/2 — ff which achieves a fidelity Fi = i(l -t- 
l/-\/2 — 77^). The purification succeeds with probability 
Pi = V^- 

For N — 2 the maximum fidelity F2 is given by 



3-772 



(48) 



One option how to achieve F2 is to choose fc = and 
Q!2,o = 77/^3 — 77^. The second alternative is fc = 1 and 
a2,i = 771/3 - 7^2/(2 - ?7^). The Hilbert space K. of the 
admissible optimal operations E is thus two-dimensional 
and spanned by basis states 



|ei) = |00)|0)-l-a2.ol*+>|l>, 
|e2> = |*+>|0)+a2,i|ll)|l). 



(49) 



To find E that maximizes P we must solve (|17|l . 

It follows from (jlTTl and that the optimal E is di- 
agonal in basis |ei), [62), E — piCi +P2e2- Consequently, 
the semidefinite program H17|) reduces to a linear program 
and we have to maximize 



P 



1 4,1 2/0 2\ 

-jPlV + -jP2il (3 - ?? ) 



under the constraints 



0<pi < 1, 



<P2 < "2,11 



< Pia2o +P2 < 1- 



The optimal aN.ki'u) that maximizes the fidelity of the 
purified state H46|l with respect to the original pure state 

(|0>+e'^^|l))/V2isgivenbyajv,fc = V4^V4+U+i and 
the fidelity of purified qubit reads 



F, 



N,k 




'k,k+l 



AN)(N) 



(47) 



The maximum achievable fidelity can be found as a max- 
imum over all fc, -FV,max = max^ -Fjv.fc- Based on numer- 
ical calculations we conjecture that for odd N the best 
fidelity is reached for fc = (iV — l)/2 while for even N 
there are two alternatives leading to the same optimal 
F, namely fc = N/2 - 1 and fc = N/2. For iV < 10 we 
have checked that these choices of fc are optimal which 
supports this conjecture. 

We shall now present explicit results for N = 1,2,3. 
Besides of maximal fidelity, we are interested also in the 



For 77 < 77th = (7 — ^/l7y/^/2 the optimal coefficients 
read pi ~ 0, p2 — 1 while for 77 > 77th we have pi — 
(•^2,1 ~ 1)/(Q;2,oQ^2,i) and p2 — a^i. For aU < 77 < 1 
the optimal probability is given by a simple formula 



P2 = lvHi~V^)- 



(50) 



Finally, when Alice sends three qubits to Bob {N — 3), 
then the optimal fidelity of Bob's purified qubit is given 
by 



5-27/2 



(4-772)^2-7/2 



(51) 



and the only way to reach F3 is to choose fc = 1 and 
aa^i = 77/1/2 — 772. The purification succeeds with prob- 
ability P3 = 77"^ - 77^/4. 

The dependence of the optimal fidelities on 7/ is plotted 
in Fig. 2(a) which clearly illustrates that the purification 
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FIG. 2: (a) Maximal fidelity of the purified qubit when Alice 
sends copies of the qubit via amplitude damping channel 
parametrized by r\. For comparison, the curve labeled -Fo 
displays the fidelity of the single qubit after passing through 
the channel, _Fo = (1 + J7)/2. (b) The corresponding maximal 
probability of successful purification. 



results in a significant improvement of the fidelity. The 
relative improvement is maximal when 77 — > but this 
is reached at the expense of very low probability of suc- 
cess, see Fig. 2(b). Note also that lim^^i P3 = 3/4. If 
Bob possesses three noisy qubits and tries to extract one 
qubit, then his optimal probabilistic strategy will have a 
finite probability of failure for arbitrarily low damping. 



V. CONCLUSIONS 

In this paper we have investigated the optimal prob- 
abilistic realizations of several important quantum- 
information-processing tasks such as the optimal cloning 
of quantum states and purification of mixed quantum 
states. We have derived a simple formula for the maxi- 
mum achievable average fidelity and we have provided an 
explicit prescription how to construct a trace-decreasing 
CP map that reaches the fidelity F^^x- We have demon- 
strated that the fidelity of probabilistic cloning can be 
strictly higher than the maximal fidelity of deterministic 
cloning even if the set of the cloned states is linearly de- 



pendent and continuous. However, it should be stressed 
that this improvement in fidelity is achieved at the ex- 
pense of a certain fraction of unsuccessful events when 
the probabilistic transformation fails and does not pro- 
duce any output state. 

The optimal probabilistic maps may find a variety of 
applications. For instance, the phase covariant cloning 
is an efficient attack on several quantum key distribu- 
tion protocols. In particular, the 2 — s- 3 phase-covariant 
cloning is explored for eavesdropping purposes in Ref. 
[2^ 1. Thus, the probabilistic phase-covariant cloning dis- 
cussed in the present paper may be possibly used as a 
new eavesdropping attack. Moreover, the general theory 
of optimal probabilistic transformations developed in the 
present paper has much broader range of applications 
than just cloning. In particular, it provides a method 
to engineer optimal protocols for purification of mixed 
quantum states. 

We have seen on the example of the amplitude damping 
channel that the optimal probabilistic purification may 
result in a dramatic improvement of the fidelity of the 
final Bob's state with respect to the original state that 
was sent to him by Alice via a noisy channel. However, 
the large improvement of the fidelity is typically accom- 
panied by a very low probability of success. It is therefore 
highly desirable to optimize the probabilistic transforma- 
tion also with respect to the average success probability 
which leads to a semidefinite program that can be very 
efficiently solved numerically. For the particular cases of 
purification of mixed states investigated in the present 
paper, we have been able to solve the resulting SDP an- 
alytically, by exploiting the symmetries inherent to the 
problem. 

The protocol considered in the present paper can be 
even further generalized as follows. One can imagine a 
scenario where the average fidelity of the operation F 
is maximized for a fixed chosen average probability of 
success P, or vice versa, these two alternatives are clearly 
equivalent. Generally, there will always be a trade-off 
between P and F and the optimal fidelity will be some 
function of P. One can then choose the working point on 
the F{P) curve that is most fitting for the particular task 
at hand. The determination of maximal F obtainable for 
some fixed P can be formulated as a semidefinite program 
similar to that given by Eq. (|17|l . The deterministic 
machines and the probabilistic machines that achieve the 
maximum possible fidelity represent two extreme regimes 
of such a more general scenario. 
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